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ABSTRACT 

A theoretical prediction to determine the frequency response of 
a cascaded remote pressure transducer system has been developed. 

This prediction was used to find a theoretical configuration which would 
produce a fairly flat frequency response (1.00 i 0.05) over a 100 hertz 
range. Such a configuration would allow cross -correlation of an unsteady, 
random, low-speed flow field by simple analog methods. 

The theoretical prediction did not hold true for configurations of 
three serially connected tubes, but the parameters, namely, transducer 
volume, tube length, and tube diameter, proved the same trends experi- 
mentally that were predicted theoretically. By observing these trends, 
a solution of 1 . 00 ± 0 . 09 over a 95 hertz range was obtained . 
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S Strouhal number = — 

T* Total temperature inside tube 

T_ Mean temperature inside tube 

o 

t Time 

T Amplitude of temperature disturbance inside tube 

U Free stream velocity 

u Complex velocity component in axial direction inside tube 

u Amplitude of U 

V Complex velocity component in radial direction inside tube 

v Amplitude of v 

V Volume of pressure transducer 
v 

V Tube volume 

Greek Symbols 

OL 
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A 
A 

CO 
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p- 
p 

T 



Wave number 

% 

Ratio of specific 
Phase angle 

Coefficient of thermal conductivity 
Absolute fluid viscosity 
Circular frequency, radians per second 
Fluid density 
Mean fluid density 

Amplitude of fluid density disturbance 
Time constant 
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I. INTRODUCTION 



Because of the increasing interest in measuring oscillatory and 
random unsteady flow fields on lightweight aeroelastic models, a need 
has gradually developed over the past several years for a direct measuring 
remote pressure transducer system to allow rapid and economical deter- 
mination of auto- and cross -correlation type information plus the attend- 
ant spectral distribution. Such a system has been developed by Bergh 
£lj at the National Aeronautical and Astronautical Research Institute 
(N.L.R.) in Amsterdam, Holland, utilizing precise lengths of tubing to 
connect the remote pressure transducers with static orifices located on 
the model surfaces. Johnson has confirmed Bergh' s one tube theory 
and Allen ^3^ utilized a two tube configuration: both of these systems 

required compensation for the system's transfer function to be included 

* > 

in the data reduction process. 

Should a configuration be found which would produce a flat fre- 
quency response over a specified range, cross -correlation of an unsteady, 
random signal could be achieved directly through the use of an analog 
circuit. Since a direct analog process does not readily allow for cor- 
rection of the system's transfer function, the primary purpose of the 
subject program was to attain a useable flat frequency response of 
1.0 — 0.05 over a 0 - 100 hertz frequency range. This would allow 
direct evaluation of the unsteady aerodynamics in many low-speed flow 
fields of interest. 
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II. EXPERIMENTAL PROCEDURE 



A. GENERAL 

All experimentation was performed either in the Aero-Mechanics 
Laboratory, Room 022, Halligan Hall, or with the IBM 360-67 digital 
computer available at the U. S. Naval Postgraduate School. 



B. THEORETICAL BACKGROUND 



Bergh [l] has developed a theoretical prediction for the pressure 
ratio of not only a single tube remote pressure transducer system 
(Figure la) but also for a cascaded tube system employing N tubes 
such as shown in Figure lb. The solution is based on several assump- 
tions; the perturbations are of a sinusoidal nature imposed on the average 
values of the parameters, fluid flow in the tubes is assumed laminar, 
and the tube lengths arq much larger than their diameters. By introducing 
these assumptions, the governing equations, i.e. , the Navier-Stokes , 
continuity, energy, and state equations are reduced to a set of partial 
differential equations which may be solved subject to satisfying the 
boundary conditions. The mass flow rates of the air in the system 
between any two adjacent tubes must be equivalent. This relation pro- 
vides an equation which reduces to the following form: 
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(b) MULTIPLE COMPONENT PRESSURE TRANSDUCER SYSTEM 



FIGURE 1 

SINGLE AND CASCADED PRESSURE 
TRANSDUCER SYSTEM 
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where: 




When a single tube configuration (see Figure la) is used, the entire 
third term of equation (1) disappears, simplifying the equation to the 
form below: 



The subscript j refers to the tubes being analyzed in the N tube system. 
From equation (1) it is seen that the pressure ratio prediction is a func- 
tion of tubes j-1 and j + 1 as well as tube j, which in turn presents diffi- 
culties when optimizing the problem for the desired result. The solution 

* 

complexity is not obvious at first; however, a closer inspection bears 
out that all the Bessel functions are functions of the imaginary parameter, 

Q /O 

i which are sometimes called Kelvin functions [4^ and hence, the 
solution is a complex quantity. The complex pressure ratio specifies 
not only the pressure ratio magnitude but also the phase shift. The 
equation also contains two uncertain quantities; the increase in trans- 
ducer volume due to diaphragm deflection, and the polytropic gas 
constant, k. Bergh proved both of these terms to have only minor effects 
on the solutions to equations (1) and (2) and hence adiabatic flow was 
assumed while transducer volume changes were ignored. The system 




( 2 ) 
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can be represented by the acoustic analogy provided by Olson CO for 
a multi-degree of freedom system. The analogy develops equations of 
motion for a sound wave of pressure p travelling in a tube and impinging 
on' a cavity volume, with the resistance represented by the viscosity and 
the capacitance term dependent on the cavity volume. Similarly, Bergh's 
solution predicts the complex pressure ratio with the resistance and 
capacitance again represented by fluid viscosity and cavity volume, 
respectively . 

C. COMPUTER TECHNIQUE 

With the aid of Bergh's equations, the pressure ratio at a certain 
frequency can now be predicted. However, because the pressure ratio 
of the next higher tube is part of the predicted ratio and the tube size 
and cavity volume interact in an irregular manner, the total solution 

does not lend itself for routine optimizing techniques to produce a flat 

* 

frequency response. Therefore the equations were programmed for a 

\ 

digital computer to facilitate pressure response prediction using pro- 
gram-operator chosen configurations. The program was designed for 
use on the IBM timesharing system but can also be used for batch pro- 
cessing, and contained a program option allowing the selection of 
either a one-, two-, or three-tube configuration. Each tube stage in a 
multiple tube configuration is specified by three parameters: transducer 
volume, tube length, and tube diameter. The program was verified for 
proper operation by comparing it with both a hand-worked calculation 
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for a one-tube configuration and Bergh's £lj results for both one- and 
two-tube configurations. 

D. EQUIPMENT 

1 . Pressure Transducers 

The two pressure transducers were designed by 
Professor L. V. Schmidt [6] , and consisted of a Bentley Detector 
System mounted behind a 0.003-inch thick annealed brass diaphragm 
which was inserted in an aluminum housing. The diaphragm thickness 
was selected to give optimum sensitivity in a 1 0.75 psi range. The 
detector system consisted of a Bentley 316 AJW proximity detector 
inserted into a brass diaphragm holder. The proximity detector was a 
reluctance gage which sensed change in the air gap between the reluct- 
ance coil and a metal target, in this case the brass diaphragm. The 
diaphragm deflection was directly proportional to an applied pressure 
perturbation, and the proximity detector output was linear with respect 
to the diaphragm deflection; therefore, the output was a relative measure- 
ment of the pressure applied to the diaphragm. 

2 . Calibration Chamber 

The calibration chamber was the same as used by Johnson 
DO and consisted of an aluminum block with a 1.5-inch diameter 
cylindrical chamber 6 inches long. The chamber ends were designed to 
accept either the transducer holders or the adapter plates which connected 
the chamber to a remote transducer through the transmitting tubes. The 
chamber had a pressure port to allow the mean static pressure to be 
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varied from an external source, and also a threaded orifice designed to 
accept either a solid plug or a University Sound ID-75 acoustical driver 
unit. Both orifices were located at the chamber's midpoint so that 
unsteady pulses were introduced symmetrically to both transducer 
systems, namely the reference transducer at one end and the complete 
transducer system at the other end. 

E. CALIBRATION 

1 . General 

All calibration procedures were identical to those performed 
by Johnson w and Allen [3], 

2 . Proximeter Calibration 

The Bentley 316 AJW proximity detectors were usually cali- 
brated in order to determine the voltage output with respect to the air 
gap between the reluctance coil face and a metal target. By mounting 
a metal target on a micrometer stem, the target could be backed off 
known distances from the reluctance coil and the output monitored on 
a D.C. voltmeter. The Bentley D252 carrier amplifiers have a variable 
amplification capability which may be adjusted until the proximity 
detector output varies proportionately with the target distance. Because 
of the large apparent sensitivity of the proximity detectors, this pro- 
cedure was foregone and the detectors were calibrated with respect to 
the pressure differential across the diaphragm during the static calibra- 
tion procedures. The proximity detectors were inserted into the brass 
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diaphragm holders until their output was -3.0 volts D.C.: this was the 



range where transducer output and diaphragm deflection varied linearly. 

3 . Static Calibration 

The transducer output was connected to a Bentley D252 
carrier amplifier powered by a -18.0 volt D.C. regulated power supply. 
By inputing the carrier amplifier output simultaneously with a biasing 
voltage into an Astro-Data D885 differential amplifier, the differential 
amplifier output, monitored on a D.C. digital voltmeter, was adjusted 
to read zero for zero pressure differential across the diaphragm. A 
U-tube manometer and hand crank pump conveniently allowed known 
pressure differentials to be introduced into the calibration chamber for 
static sensitivity determination. The complete instrumentation diagram 
is shown in Figure 2. Each proximity detector was calibrated during 
this stage by introducing positive and negative pressure differentials 
of equal magnitudes into the calibration chamber and adjusting each 
carrier amplifier until the output was approximately the same for equal 
but opposite pressure differentials over a -0.8 to +0.8 psi range. The 
static calibration curve is shown by Figure 3 and bears out the linear 
sensitivity variation as measured at the carrier amplifier output of 1.37 
and 1.04 volts/psi for transducers 1 and 2, respectively, over a ±0.8 
psig pressure range. 

4 . Dynamic Calibration 

For the dynamic calibration, a sinusoidal perturbation was 
introduced into the chamber by closing the static pressure port and 
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FIGURE 2 

STATIC CALIBRATION INSTRUMENTATION 
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STATIC CALIBRATION CURVE 
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inserting an ID 75 University Sound acoustical driver unit driven by a 
Wavetek 114 signal generator which was amplified by a power amplifier. 
The acoustic driver unit was adjusted to operate at an input level of 
13.5V RMS which corresponds to an oscillating pressure amplitude in 
the calibration chamber cavity of approximately 0.13 psi at 120 hertz. 
The carrier amplifier output was monitored with both a true RMS meter 
and a precision phasemeter, and visually monitored on a dual beam 
oscilloscope. The dynamic calibration instrumentation diagram is shown 
in Figure 4. The frequency was varied from 10 - 120 hertz in 5-hertz 
increments. Since the disturbance was introduced symmetrically into 
the chamber, there was no appreciable phase shift in oscillating 
measures at either end of the chamber; however, the sensitivity of 
transducer #2 was 0.787 of transducer #1 at every frequency increment. 

F. EXPERIMENTATION 

The data runs were made with the same system as was used for 
the dynamic calibration procedure with one exception; transducer #2 
was remotely located from the calibration chamber through the use of 
nylon and stainless steel tubing. The brass diaphragm had a centered 
circular wrinkle which allowed brass shims to be inserted between the 
diaphragm and the plexiglass spacer in the transducer housing without 
hindering movement of the diaphragm. The wrinkle volume was estimated 
as 0.00211 in. and the shims had an inside diameter of 0.523 in.: 
this allowed any desired transducer volume to be constructed. Frequen- 
cies of 10 - 120 hertz were investigated in 5 -hertz increments for 
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different one-, two-, and three-tube configurations. The ratio of trans- 
ducer output #2 divided by transducer output #1 was calculated and, 
because the sensitivity of transducer #2 was determined during dynamic 
calibration of transducer #1 as 0.787, the quotient was divided by 
0.787 to normalize the output of transducer #2 with respect to transducer 
#1. This determined the frequency response (transfer function) of a given 
configuration . 
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III. DISCUSSION AND RESULTS 



The computerized solution of Bergh's [i] prediction allowed 
quick and easy determination of the frequency response for a one-, two- 
or three-tube system, and could be easily modified for a greater number 
of tubes . 

The one-tube solution generally followed the typical organ pipe 
solutions which produced sharp resonance peaks at many frequencies. 
When the resonance peaks did not occur, a heavily damped solution 
appeared. A solution was sought that would reduce the resonance peaks 
or raise the damped solution to the desired goal (1.00 ± 0.05); however, 
no such solution could be found for a 0 - 100 hertz range. The data runs 
did confirm Bergh's theoretical predictions for both pressure ratio and 

phase shift for a one-tube configuration, as can be seen by Figure 5. 

** 

Therefore, a two-tube solution was examined in hopes that the pressure 
of an additional tube would counter the undesirable characteristics of 
the first tube. Such a configuration was found with the desired frequency 
response of 1 . 00 ± 0 . 05 from 10 to 100 hertz as shown by Figure 6 . 
However, it was felt that the first tube in the system had a diameter 
that would result in the static orifice on the body being large enough to 
perturb the flow in the region of static pressure measurement. The two- 
tube configuration also verified Bergh's theoretical prediction of phase 
shift, but the pressure ratio was consistently lower (6%) than predicted 
for frequencies greater than 40 hertz (see Figure 7). A three-tube 
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PHASE SHIFT, DEGREES PRESSURE RATIO, P 2 /P 




FREQUENCY, HERTZ 

0 20 40 60 80 100 120 




TRANSDUCER VOLUME (CU. IN.) 0.0195 
TUBE LENGTH ( IN. ) 28.0 

TUBE DIAMETER (IN.) 0.078 

FIGURE 5 

SINGLE TUBE CONFIGURATION 
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DESIRED THEORETICAL SOLUTIONS 
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FIGURE 7 

TWO TUBE CONFIGURATION 
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configuration was then sought in order to reduce the diameter of the first 
tube to an acceptable value (less than 0.05 inch). A solution which met 
all desired specifications (1 . 00 ± 0 . 03 over a 10 - 100 hertz range) was 
found (see Figure 61 and physically constructed in order to obtain experi- 
mental results to confirm the configuration selected. The configuration 



specifications were: 


Tube 


Tube 


Tube 




# 1 


# 2 


#3 


Transducer volume (in. 3) 


0.0 


0.0 


0.0238 


Tube length (in.) 


0.45 


29.0 


21.0 


Tube diameter (in.) 


0.0465 


0.078 


0.03 



The three-tube solutions again confirmed the phase shift prediction; 
the experimental values were within 7 % of their predicted values (see 
Figures 8a-d). However, as was the case with the two-tube solution, 
the three-tube experimental configuration produced pressure ratios less 
than those predicted by Bergh's solution for frequencies above 30 hertz 
with the values lower by as much as 30% (see Figures 8a-d). This 
characteristic made the theoretically acceptable solution (Figure 8a) of 
limited practical value since the experimental response was consider- 
ably below the prediction. Even though the actual solutions were not 
as predicted, Figure 10 shows trends comparable to those predicted by 
Bergh's solutions for identical configurations (Figure 9). 

From Figures 8b-d, it was observed that changing the lengths of 
tubing raised or lowered the overall curve and a decrease in transducer 
volume (compare Figures 8a, b) tended to increase or raise the overall 
curve in the region of interest and decrease the slope above 80 hertz. 
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PHASE SHIFT, DEGREES PRESSURE RATIO 
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0.0 0.0 0.0238 
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0.0465 0.078 0.03 



FIGURE 8a 

THREE TUBE CONFIGURATION 
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PHASE SHIFT, DEGREES PRESSURE RATIO, P 4 /Pj_ 
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TUBE LENGTH (IN.) 
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FIGURE 8b 

THREE TUBE CONFIGURATION 
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PHASE SHIFT, DEGREES . PRESSURE RATIO, P 4 /P 
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FIGURE 8c 



THREE TUBE CONFIGURATION 
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FIGURE 8d 

THREE TUBE CONFIGURATION 
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TRANS. VOL. (IN . J ) 0.0238 0.0238 0.0174 0.0174 

# 2 LENGTH (IN.) 31.0 29.0 29.0 28.5 

# 3 LENGTH (IN.) 21.0 21.0 21.0 20.5 
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THEORETICAL TRENDS 



TRANS. VOL. (IN. J ) 0.0238 0.0238 0.0174 0.0174 

# 2 LENGTH (IN.) 31.0 ■' 29.0 29.0 28.5 

# 3 LENGTH (IN.). 21.0 21.0 21.0 20.5 
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EXPERIMENTAL TRENDS 



A small volume was tested and a possible solution of 1.00 ± 0.09 over 
a 10-95 hertz range was finally attained (see Figure 8d) . This configu- 
ration appears acceptable for direct analog data processing in low-speed 
fields of interest. Much of the separated flow phenomenon is involved 
at non-dimensional frequencies or Strouhal numbers, S, less than 0.4 
which would typically correspond to frequencies less than 100 hertz for 
a five-inch diameter circular cylinder at a tunnel speed of 100 feet per 
second. The solution configuration is specified below: 





Tube 


Tube 


Tube 




#1 


# 2 


# 3 


Transducer volume (in. 3) 


0.0 


0.0 


0.0174 


Tube length (in.) 


0.45 


28.5 


21.0 


Tube diameter (in.) 


0.0465 


0.078 


0.03 



In addition to the tube lengths and transducer volumes, the three 

tube diameters were very important parameters when attaining the theo- 

* 

retical response desired. The diameter effects and trends were not 
examined, however, due to time and facility limitations and particularly 
due to the limited selection of commercial tube diameters available with 
acceptable inside diameter specifications. 

The purpose for attaining a flat frequency response was to allow 
easy auto- and cross -correlation of a body in a low-speed unsteady 
flow employing a simple, low cost analog method shown by Figure 11. 
The cross-correlation coefficients, which relate the coherence of two 
physical measurements, say at stations (1) and (2), are in general 
defined for a specific value of time lag, T , by: 



35 



+15 1 MEG 



o 

o 

o 




■O CQ 
+ 




o 



o- 



■o 



-i 

< 

CX 

U4 

H 

X 

UJ 



o 



>- 

i 

X 







o 



<D 



CNJ 

<D 



o 

o 



w 

( X 

O 

K— < 

Ph 



36 



CROSS-CORRELATION CIRCUIT 



n /o^\ _ (J) ^2. 

However, the elementary analog correlation circuit shown in 
Figure 11 does not have a time lag capability and hence the benefit of 
simplicity is associated with a loss of information. Therefore, the 
information "to be obtained must be at zero time lag as specified by 
equation (3): 

p ,VI_ <e, (°) &-U M > 

K,lU <e, l >''2 <e£> v i 

Since the system in consideration, i.e. , a remote pressure trans- 
ducer system, can be likened to a typical transfer function system shown 
by Figure la , b , the cross-correlation inputs, e^ and e^, can be repre- 
sented by the transducer outputs, p^ and p^. If the disturbance were 
of a sinusoidal nature as it was for the laboratory experiment, then the 
two measuring stations become strictly functions of the frequency. 

Suppose: 

7°i - 7°oi sin co t 

7^- 7°^ Sw (cot- 0) 

then: 

■Pi (+)= 1°o\ &, LoS) sin (co +-<£,) 

-^ i (+)= fw (^) Svn (cut- ©- 
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Substitution into equation (3) will yield a cross -correlation coefficient 



as follows: 



R 



12 . 



(o') - — ^ ^ ^ ^ 



G t~ (<*■>) 



From equation (4) it can be seen that for the case of a pair of information 
signals at a unique circular frequency, co> with the respective transfer 
functions {co) and (c*>) equivalent, then: 

R l% (°) “ CC>S & 

and a flat frequency response is not necessary. It is only necessary to 
have identical transfer functions used with both signals. 

However, this analysis eliminates one of the most important types 
of flow fields, the unsteady random signal such as that produced by flow 
separation from a bluff body. In such cases the cross -correlation 
coefficient can no longer be simplified to equation (5). Instead, because 
of its reliance on G as a function of CO , the output of the cross- 
correlation system becomes: 

£ £ v _ JT 

f [ X §, [ JT 3 4 (6) 



where 




power spectral density of signal 1 

power spectral density of signal 2 
co-spectral density of signals 1 and 2 
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At the input of the transducer system, the cross -correlation coefficient 



is not a function of G and thus reduces to the form: 



r 12 .(°)= 



JT 






(7) 



In order for (6) to equal (7) for any random signal, the following relation- 
ship is obviously necessary. 



G^ - Gg (ct>) - l'O 



( 8 ) 



Should (8) be the case, then, 



R 12 (0) input R 12 (0) output 



and the system no longer relies on the gain of the transducer system. 

A complete derivation of the cross -correlation equations may be found 
in Appendix C. 

Because of the flat frequency response, the collected data may be 
inputted directly to the cross -correlation circuit shown in Figure 11 and 
the output monitored by either a D.C. voltmeter or an x-y plotter to 

i i 

2 2 2 2 

obtain quantities ^e^e^ ' ^ e l ^ ' anc * ^ e 2 ^ a ^ ows a simple 

calculation of the cross-correlation coefficient at zero time lag. 
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IV. CONCLUSION AND RECOMMENDATIONS 



Although the breakdown of Bergh's theory seems to be directly 
proportional to the number of tubes in the configuration, it does predict 
fairly accurately the trends that each parameter impresses on a system. 
There is, therefore, a good possibility of obtaining a flatter frequency 
response than the one obtained in Figure 8d . Time limitations required 
the termination of the experiment and the three-tube configuration was in 
no way completely examined. 

It should be noted that the flat frequency response traits considered 
herein applied to a fairly restrictive range of frequencies, namely less 
than 100 hertz. The processing of a broader class of signals, i.e., a 
higher range of frequencies, will probably require the utilization of several 

steps or procedures, including: 

* 

a. Recording of the analog form of the data signals on a magnetic 
tape recorder. 

b. Digitizing the analog signals using a common time base or 
clock, and 

c. Processing the digitized data by means of a digital computer. 

Portions of this procedure have been within the scope of the school's 

experimental facilities within recent times, cf. , Jesberg ^7^ ; however, 
new equipment additions or replacements make this former technique no 
longer practical until a repeat is done of the work by Jesberg ^7} . 
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It is observed from Figures 7 and 8a-d that the theory under- 



predicts the pressure ratio for both two- and three-tube configurations, 
yet the phase prediction is very close to the correct value. 

The following recommendations are therefore made: 

1. ) Further experimentation of parameter effects in order to 

obtain a flatter frequency response. 

2. ) Possible modification of the computer program to allow 

compensation for the difference between the actual and 
theoretical response. 
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APPENDIX A 



THEORETICAL RESPONSE OF A 

REMOTELY LOCATED PRESSURE TRANSDUCER SYSTEM 

Bergh's theoretical model applies to any number of transmitting 
tubes in a cascaded tube system as shown by Figure lb; however, many 
assumptions are made in order to reduce the governing equations to a 
workable form. These assumptions are as follows: 

a) the sinusoidal disturbances are very small: 

b) the inner diameter of the tubing is small with 
respect to its length: 

c) laminar flow exists throughout the system: 

d) the unsteady disturbances are perturbations 
imposed on the average values on the parameters 
indicated below: 



1) 


P = 


P o 


+ 


iwt 

pe 


Pressure 


2) 


P = 


P o 




iwt 

pe 


Density 


3) 


T = 


T 

o 




m iwt 
Te 


Temperature 


4) 


U = 


u 

o 


+ 


iwt 

ue 


Axial Velocity 


5) 


V = 


V 

o 


+ 


iwt 

ve 


Radial Velocity 



These assumptions allow the governing equations to be reduced to the 
following forms: 
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Navier-Stokes Equations 



_L 




gfru 



L ^ u - 

Y 9 r 




A co u. 




Continuity Equation 




Perfect Gas Law 




Energy Equation 

^ [ |r- + 1 U .] r + ^ «<«».% c j* T 



Since the equations are now partial differential equations, 
conditions can now be specified: 

a) At r = R, the tube's inside wall, 

1) u(X, R) = 0 

2) v (X , R) = 0 

3) T (X , R) = 0 



(A2) 

(A2a) 

(A3) 

(A4) 

(A5) 

boundary 
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b) At r = 0 , the tube center, 



1) v(X, 0) = 0 

2) u(X,0) = finite 

3) T(X,0) = finite 

4) p(X,0) = finite 

5) p(X,0) = finite 

From equation (A2a) it is evident that p = p(x) , which allows a 
separation of variables technique. 



assume: 

T= foo K(a) 1 

where: 

, ( PoD\V Z. 

0(3 ^ ) 

p K - - PRflWOT^ M umber 

A . 

2 = T p ' 

Equation (A5) takes on a Bessel equation form of the following type in 
this type of solution: 



h (b) = c, S 0 (a) f c*. Y c (a) +■ 



-p / 

f(*) 



(A7) 



Apply boundary conditions a3 and b3 



fOO* 



i -pOO 

c, f>o%c P 



(A 8) 
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Using (A7) and (A8) in (A6) yields: 




O'. («. P,'^) J j 



(A9) 



Similar methods using boundary conditions al and b2 can be applied to 
solve equation (1) as given below: 



U- 




i o;(«<) 




(A10) 



Substituting equations (A91 and (A10), integrating with respect to r, 
and applying boundary condition a2 will put equation (A 7) in the form: 




S’, . v^t ay (<* v*" 1 ) 



R> <J Z P Sz.(A 

2- dx 1 (=<•) 



Since v(X,0) = 0 implies F(X) = 0, 



60*- 

a? 




3k 

J " <«. A- ‘^) 






JL (°0 c/ 2 -^ __ _ 
Jo (<*•) c/x^ 



The above equation can be solved for the homogeneous solution: 

-p= + B='* x 



where 

j CO 1 * Jo 




I 

i. 



4* 






Ji (<*• P^) 

Jo (<* fW 
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The coefficients A and B can be determined only after the boundary 

conditions for both ends of the tube have been specified. It should 

also be noted here that the Bessel functions in the equations are complex 

3/o 

numbers since c*. = c< (i y). This is discussed further in Appendix B. 

Two additional assumptions must now be made if the solution is to 
cover an N-tube, N volume configuration: 

c) ’ the pressure and density in the instrument volumes 

are time dependent only 

d) the pressure expansion in the instrument volume 
is a polytropic process described by 

■py 

— = constant 

K 

The constants A and B are now determined in the following procedure. 

* 

For any tube j , 
a t x. = 0 

J . (All) 

p. = A. + B. 

J-l J J 




at x. = 1 . 
J J 



-Pi - 



n. - iii 
3 



Jo ( 

Jo OH) 



-I 



B e -A 

0 o 



“1 



(A12) 
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Equations (All and (A12) allow solution for A. and B., 

8 . = Pj-i 

* - e~ 



The mass leaving the tube can now be calculated: 

/M >, = p H u i ^try-dr = 



IT** fcJlC**) r 



CO ) i _ ^ 






For any tube j + 1 , 



at x. . = 0 
J+l 



p. = A. , + B. , 
J J X 1 J-*-l 



at = L 



>+' 



■Pj+i = A*+t f B , +1 e -^ |L > +l 



J + L 






] 



The mass leaving tube j+l is now calculated: 

^#0= /«.***' ^ “i. ^irrdy- = 



ILiill ^ a, '~ 1 ' 1 J i- (: °<vn') (Aj + |- 8 i+ , ) 

■ 3 ” 5 C 0 ^**) 



J. Uj 



(A13) 



(A14) 
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Equations (A13) and (A14) allow solution for and 

_P£+I -piC~' t ’ i *' Li+ ' 



g4i+ii£+i ^ 



B 




- -pin 

(ofo+iLs+t _ (o-tiul-i+t 



The mass variation within the instrument volume must equal the differ- 
ence in the mass leaving tube j and entering j+1, namely: 



d VYlv -** tu y 

Jt = C ft *i 





(A15) 



Substitution of the pressure coefficients into (A15) and recalling the 



definitions 



cosh cot = 






and sinh CO t 






the final pressure ratio prediction is obtained. 

tl 



ft ,- H 



-f ILEjdti ‘fta-H 310**1) Jo s/nh x 

^ ^ <K -JoC^+i) J 16 *$) s >m ( ta+, L i+') 

| COSW + "* 
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APPENDIX B 



COMPUTERIZATION OF THE THEORETICAL RESPONSE 
OF A REMOTE PRESSURE TRANSDUCER SYSTEM 



A. Computer Program Development 

In Bergh's theoretical development (Appendix A) , the Bessel 

3/2 

function parameter, , contains the imaginary term , i , which can 
easily be represented by the Kelvin functions Ber, Bei, Ber^, and Bei^. 
These functions were used in the program in the form of an infinite 
power series and terminated when the terms become smaller than 10 
The series expansions were taken from McLachlan £f] and altered into 
the form below: 



vTo6*) - 
cX - 




(Bl) 



be \r y 




Y+ , y* 

(Z'f) 7 ' 4^6-e) 2 - 



V 



a 



( Z"f'6 • 8. to >12.)' 



•K.. 



be i 




y C y l o 

■ . — l 4- 1 . . 

(~L (1. .• 4' t 



49 



0*) = ber, y + J- be l , y 




(B2) 




/ r v 1 y* . y 1 



y 



(B3) 



B . Computer Program Operation 

The computer program was designed primarily for use on a time- 
sharing system remote terminal where quick output could be obtained 
for a given input. However, it may also be modified for batch processing 
as explained below. The computer program allows selection of either 
a 1, 2 or 3 tube configuration by the operator and prints out the pressure 
ratio and phase shift for multiples of 5 hertz from 10 to 120 hertz. The 
flowchart is given in Figure (Bl). 
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START 





FIGURE B1 

COMPUTER PROGRAM FLOWCHART 
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APPENDIX C 



CROSS-CORRELATION OF A RANDOM UNSTEADY SIGNAL 



The cross -correlation coefficient may be defined in its most 
general form similar to methods described in Bendat and Piersol DO. 

r,Y(+) - <£'<♦> 

However, because of the zero time lag limitation of the analog circuit 
shown in Figure 11, this reduces to the form: 



R,z(°) 






(C 2) 



For the system input, the following terms are defined as: 

e, (+) e z (+) d-t- 

= f<T 

w>= c £ u Mc/co 

^ = i'c (‘‘-O «/co 



(C3) 
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in 

i 22 



= power spectral density of signal 1 
= power spectral density of signal 2 
= co spectral density of signals 1 and 2 



After passing through the system, the quantities have a transfer function 
imposed on their true output, and the following terms are defined: 



e'.M = 6,(0.) e e, (o.) 
e,M = A,e-^° + 



(C4) 



thu s : 

> = _y o °°6,(o.) C^) dco 

5o (^(o.) jf,, (.L**) d to 

= £» c£<««) «)«/« 



A relationship between the output and input cross-correlations may be 



obtained by direct substitution of (C3) and (C5) into (C21 , yielding: 

/\ 

• oo 



Jo 5,^Ct>)c/to 



^ J J z c<*>) <j uTJ 

jT &,<») jjT |t (o>) c/co 

{pr f „ (-) ] [5 ” If h ^ ^] ? ,/i 



(C6) 
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There are two possible simplifications available at this point. If the 
input were of a harmonic nature, say ac frequency Co , then the expres- 
sions for spectral density would simplify and yield: 

§ (uj-Uj,) J CO ■= 



(C7) 



similarly: 

<.§ 1 2 >=f,r sctwjjtu =i„cw,) 

<e i 2 > := jT" c< o) 

(C8) 

<£, £,> = 5Tl,^) <V“) 

and the dependence of the zero time lag cross -correlation coefficient 
upon the system's transfer function would disappear due to cancellation. 

However, if the input were of a random nature, then the only 
simplification readily available is to have: 

G x (*>) = G 2 (co) = 1.0 

in which case the original input cross -correlation is obtained since the 
system transfer function factors out from the integrands, cancels, and 
simplifies to equations.* 

<e, l > = 5 „(<o) ^ 

< ei > - < e ^> = Sr 1^) Jlo 
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< mo *<«.*».> = dw 

Thus: 

Riz. Co) - R^-lCo) (cio) 



Since the latter case also applies to the harmonic signal data 
processing, this is a most desirable situation and is the motivation 
behind seeking a flat frequency of 1.00 ± 0.05. 
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